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ABSTRACT 

A primary goal of the FUSE niission is to understand the origin of the O VI ion in the interstellar 
medium of the Galaxy and the Magellanic Clouds. Along sightlines to OB-type stars, these interstellar 
components are usually blended with O VI stellar wind profiles, which frequently vary in shape. In order 
to assess the effects of this time-dependent blending on measurements of the interstellar O VI lines, 
we have undertaken a mini-survey of repeated observations toward OB-type stars in the Galaxy and 
the Large Magellanic Cloud. These sparse time series, which consist of 2-3 observations separated by 
intervals ranging from a few days to several months, show that wind variability occurs commonly in O VI 
(~60% of a sample of 50 stars), as indeed it does in other resonance lines. However, in the interstellar O VI 
A1032 region, the O VI A1038 wind varies only in ^--^30% of the cases. By examining cases exhibiting large 
amplitude variations, we conclude that stellar-wind variability generally introduces negligible uncertainty 
for single interstellar O VI components along Galactic lines of sight, but can result in substantial errors 
in measurements of broader components or blends of components like those typically observed toward 
stars in the Large Magellanic Cloud. Due to possible contamination by discrete absorption components 
in the stellar O VI line, stars with terminal velocities greater than or equal to the doublet separation 
(1654 kms^^) should be treated with care. 

Subject headings: line: profiles - stars: winds - stars: mass-loss - stars: early-type 



stellar components in a complicated, time-dependent way, 
which may compromise measurements of interstellar O VI 
by unknown amounts. 

To assess the influence of stellar wind variability on 
measurements of interstellar O VI, we have used FUSE 
to obtain 2-3 observations for each object in a sample of 
~50 OB-type stars in the Galaxy and Large Magellanic 
Cloud (LMC). These sparse time series exhibit variations 
over intervals ranging from about a day to several months. 
Preliminary results suggest that stellar wind variability is 
often found in O VI, with no apparent differences between 
the Galaxy and the LMC. In this Letter we illustrate the 
influence of the variability on interstellar measurements 
along five sightlines, which typify the problems that may 
be encountered. 

2. OBSERVATIONS AND ANALYSIS 

The design and performance of the FUSE spectrograph 
have been described by Moos et al. (2000) and Sahnow et 
al. (2000), respectively. Table 1 summarizes the repeated 
FUSE observations for the five sightlines discussed here. 
The observations of HD 97913 were obtained in histogram 
mode; all others were obtained in time-tag mode. All ob- 
servations were obtained through the 30" x 30" apertures. 
Standard processing with the current version of the cal- 
ibration pipeline software (calfuse vl.8.7) was used to 
extract and calibrate the spectra. The extracted spectra 
associated with separate exposures of a given observation 
were aligned by cross-correlating the positions of strong in- 
terstellar lines (other than O VI). Since these subexposures 
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1. INTRODUCTION 

The O VI AA1G32, 1038 resonance doublet is an excel- 
lent diagnostic of high-energy processes (> 113.9 eV) in 
rarefied media because of the large cosmic abundance of 
■ oxygen and the substantial energy required to strip five 
electrons from O atoms. A key program for the Far Ultra- 
violet Spectroscopic Explorer {FUSE) mission is to mea- 
sure the strength of this doublet along many sightlines, in 
order to determine the origin and evolution of hot gas in 
the interstellar medium (ISM; see, e.g.. Savage et al. 2000; 
Sembach et al. 2000). 

However, O VI also occurs in the stellar winds of the 
OB-type stars that frequently serve to illuminate the inter- 
vening ISM. Its presence in these winds is surprising, since 
the photospheric radiation fields of only the hottest stars 
produce atmospheric O VI. A current explanation is that 
winds contain ensembles of strong shocks, which might be 
produced by the strong line-driven instability (Owocki et 
al. 1988). The X-rays created by these shocks produce 
O VI from the dominant O IV ion via Auger ionization 
(Cassinelli & Olson 1979). Consequently, the shapes and 
strengths of the O VI wind profiles are largely determined 
by the distribution, strength, and time-dependent behav- 
ior of the shocks. Indeed, previous work with Copernicus 
(see, e.g., York et al. 1977; Snow et al. 1980) and early 
results from FUSE (Massa et al. 2000) confirm that the 
O VI wind lines of early-type stars are variable on a va- 
riety of time scales, as are most other wind features; see 
Prinja (1998) for a recent review. As a result, undulations 
in the broad P-Cygni wind profiles blend with the inter- 
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did not exhibit significant variations, they were coadded 
and rebinned to a spectral resolution of ^15,000. 

Fig. 1 shows the FUSE spectra of the LMC star 
Sk —69° 124 and summarizes the method used to detect 
stellar wind variability in these sparse time series. (1) 
Since the flux is well calibrated, two spectra of the same 
object taken at different times were overplotted and di- 
rectly compared (top panel). (2) To ensure that the wind 
variations are genuine, we require that they are similar in 
both transitions of the doublet (bottom panel) and in data 
from different detector segments. Since the O VI A 1032 
line is intrinsically stronger than the A1038 line, the am- 
plitude of any fluctuation will generally be larger in the 
blue component. Fig. 1 shows clearly the dramatic effect 
that the stellar wind variability might have on interstel- 
lar measurements. If the O VI A1038 Galactic feature is 
detected, the LMC component is not; but even more strik- 
ing is the absence of any obvious interstellar O VI A1032 
component. 

To quantify the influence of the stellar wind variabil- 
ity on the interstellar measurements, we consider three 
other stars that exhibit large variations in the vicinity of 
the interstellar O VI A1032 feature (Table 1 and Fig. 2). 
Typically, the signal-to-noise ratios per resolution element 
lie between 15 and 30 for these spectra. We have fitted 
the apparent stellar continuum of each observation inde- 
pendently; i.e., we did not fit the continua with a pre- 
conceived notion based on the other spectrum. For that 
reason, the chosen continua look different between pan- 
els in Fig. 2. The criterion to define the continuum was 
that the O VI A1032 full width (FW) should not be larger 
than about [—70,-1-70] kms^^ for the Galactic stars (Sav- 
age et al. 2000), while for the LMC it should be approx- 
imately in the range of [—80,-1-350] kms~^ (Howk, priv. 
comm., 2001), with small scatter from one sightline to an- 
other. The O VI A1032 line is also bracketed by two H2 
lines, (6 - 0) P(3) A1031.19 and i?(4) A1032.35, which 
often complicate continuum placement. Although a de- 
tailed H2 model based on all the transitions in the FUSE 
bandpass could in principle constrain the influence of these 
lines on the local continuum, for this study we have sim- 
ply interpolated over them. Moreover, if the interstellar 
O VI A1038 feature is resolved, it should always be com- 
pared to O VI A1032, though generally O VI A1038 is 
very difficult to discern due to blending with C II, C II*, 
and H2; this is the case for all the stars discussed here. 
The continuum was fitted by a Legendre polynomial of 
degree d (see Table 1). Fig. 2 shows very complex pro- 
files and therefore the order of the polynomial to fit the 
continuum needed to be relatively large. Once the con- 
tinuum was considered satisfactory, the equivalent width 
and column density were calculated by directly integrat- 
ing the O VI A1032 intensity and apparent column density, 
Na{v) = 3.768 X 10^^Ta{v)/[fXaiA)] (Savage & Sembach 
1991) where / = 0.132 and Aq = 1031.926. The results 
are given in Table 1. The statistical errors and errors 
on continuum placement, presented in Table 1, were ob- 
tained following the method outlined by Sembach & Sav- 
age (1992). The continuum placement errors are compara- 
ble to or larger than the statistical errors and do not reflect 
systematic differences that result from different choices for 
the wavelength intervals defining the (unknown) location 



of the true continuum. 

For the two Galactic stars, there is good agreement be- 
tween the absolute measurements of the interstellar O VI 
features, though the la errors can be quite large (see Ta- 
ble 1). However, even though the H2 ^"(3) and O VI fea- 
tures are similar for both the normalizations illustrated in 
Fig. 2, the depth of the H2 i?(4) hue is quite different, 
which suggests that the continuum is not well defined on 
the red side. The discrepancy is particularly obvious in the 
case of HD 168941. For HD 97913, the two normalized pro- 
files have noticeable differences, even though the measured 
equivalent widths agree within the uncertainties. For this 
star, we did not account properly for the CI I feature; this 
also explains why the H2 P(3) line appears to be over- 
normalized. The situation was even worse for the second 
observation and therefore we concentrated on a smaller 
part of the spectrum. The sightline toward the LMC star 
Sk —67° 191 shows large discrepancies depending on which 
normalization is adopted, so that the resultant equivalent 
widths do not agree within the estimated la uncertainties. 

3. DISCUSSION 

Stellar wind variability can have serious consequences 

for interstellar O VI studies: (1) In some instances, it 
can make the interstellar O VI identification very diffi- 
cult (see Fig. 1). (2) It complicates significantly the stel- 
lar continuum placement which increases the uncertainty 
in determining the interstellar line strength, or even gives 
erroneous results in situations where multiple interstellar 
components produce broad features. Since the interstellar 
O VI FW toward LMC stars is typically about three times 
that observed in Galactic stars, analysis of LMC sightlines 
are particularly susceptible to these uncertainties. While 
this study shows that stellar wind variability is generally 
negligible for studies of interstellar O VI along sightlines 
toward Galactic OB-type stars, systematic errors in con- 
tinuum placement can remain large for complicated con- 
tinua. Modeling of stellar wind profiles could in principle 
help to constrain the shape of the local continuum, but at 
present the fits are uncertain for Galactic stars (due mainly 
to the H2 contamination), and moreover, do not yet incor- 
porate time-dependent behavior in a self-consistent way. 
However, since these cases show wind variations very close 
to the interstellar O VI A1032 line, they should be consid- 
ered as "worst case" scenarios. Preliminary results from 
our mini-survey for O VI variability (Lehner et al. 2001, in 
prep.) indicate that at least ^60% of our sample exhibits 
O VI wind variability; but less than ~30% (14/50) of the 
total sample exhibits variations in the region of the inter- 
stellar features, i.e., —600 < f < 600 kms"'^, at the time 
of the survey. The relevant velocity ranges are defined by 
the separations between the O VI doublet and O VI A1032 
and H I Ly /?, 1654 and 1805 kms~^, respectively. For 
10 of these 14 stars, the terminal velocity (woo) is known: 
7/10 have roughly 1665 < v^o ^ 2135 kms^^, while 3/10 
have lower Vao and can blend with high velocity interstellar 
components. 

Although the nature of the stellar wind variations can- 
not be determined from the mini-survey, they arc likely 
related to the discrete absorption components (DACs) 
that are commonly observed in the wind lines of hot 
stars (Prinja 1998). DACs are particularly strong near 
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Woo- When Vryo 1650 kms ^, the interstellar mea- 
surement will not be affected by DACs (provided that 
t'iSM is not much greater than kms~^). But when 
Woo ~ 1650 kms^ , contamination by DACs is possible. 
When Woo ^ 1650 kms^^, the effect of the DACs on 
the local continuum is much less, since they change dur- 
ing their blueward motion from broad, weak features to 
a deep, sharp one near Wqo (Prinja 1998), assuming that 
O VI DACs behave as for the other ions. It is also impor- 
tant to note that the separation between O VI A1032 and 
the edge of the strongly saturated part of H I Ly can be 
<^ 1800 kms^^, and therefore information from the DACs 
in O VI A1032 can be lost. Fig. 3 shows the blueward mo- 
tion of a DAC in the wind of the Galactic star HD 91597. 
Three other sightlines presented here show an absorption 
feature on the red side of O VI A1032 (see Figs. 1 & 2). 
In the more extended sample that will appear in our later 



paper, only one case has a FW comparable to a Galac- 
tic interstellar O VI line, while all others have typically 
200 < FW < 400 kms"\ Although we are not able to 
quantify further the effects of DACs from short time series 
of many different objects, it is clear that their presence in 
O VI wind profiles can complicate continuum placement 
(see Fig. 2) and in the worst case situations even produce 
erroneous measurements. 

We thank Chris Howk and Dave Bowen for insightful 
comments on an early draft of this Letter, and Don York 
for initially suggesting a study to gauge the impact of stel- 
lar variability on measurements of interstellar O VI ab- 
sorption. This work is based on data obtained for the 
Guaranteed Time Team by the NASA-CNES-CSA FUSE 
mission operated by the Johns Hopkins University. Fi- 
nancial support to U. S. participants has been provided 
by NASA contract NAS5-32985. 



REFERENCES 



Cassinelli, J. P., & Olson, G. D. 1979, ApJ, 239, 304 

Howarth, I. D., Siebert, K. J., Hussain, G. A. J., Prinja, R. K. 1997, 

MNRAS, 284, 265 
Massa, D., et al. 2000, ApJ, 538, L47 
Moos, H. W., ct al. 2000, ApJ, 538, LI 

Owocki, S. P., Castor, J. I, & Rybicki, G. B. 1988, ApJ, 335, 914 
Prinja, R. K. 1998, in Cyclical Variability in Stellar Winds, ed. 

L. Kaper & A. W. PuUerton (Berlin: Springer), 92 
Prinja, R. K., & Crowther, P. A. 1998, MNRAS, 300, 828 



Sahnow, D. J., et al. 2000, ApJ, 538, L7 

Savage, B. D., & Sembach, K. R. 1991, ApJ, 379, 245 

Savage, B. D., et al. 2000, ApJ, 538, L27 

Sembach, K. R., Savage, & B. D. 1992, ApJ, 83, 147 

Sembach, K. R., et al. 2000, ApJ, 538, L31 

Snow, T. P., Wegner, G. A., & Kunasz, P. B. 1980, ApJ, 238, 643 
York, D. G., Vidal-Madjar, A., Laurent, A., & Bonnet, R. 1977, ApJ, 
213, L61 



4 



Lehner et al. 



Table 1 

Program star summary and interstellar O VI measurements 



Star 


Sp. Type 




UT-Date 




M 


^exp 


d" 










(kms-^) 




(days) 




(s) 




(mA) 


(dex) 


Sk -67°191 


08.5V 


1750"^ 


1999-12-17 


287.0 


2 


8460 


6 


440 ± 23 ± 53 


14.67 ±0.07 








2000-09-29 




1 


6569 


5 


609 ± 23 ± 45 


14.85 ±0.05 


Sk -69°124 


09 II 


1430'' 


2000-10-03 


62.1 


1 


2120"= 














2000-12-04 




2 


4220'^ 








HD 168941 


09.5 II-III 


1795^ 


2000-08-30 


2.2 


2 


4110 


6 


221 ± 24 ± 25 


14.46 ±0.11 








2000-09-01 




2 


3982 


7 


221 ±16 ±51 


14.49 ±0.13 


HD 97913 


B0.5IVn 




2000-05-26 


1.1 


8 


5996 


6 


195±9± 13 


14.43 ±0.05 








2000-05-27 




12 


9980 


3 


204 ±11 ±40 


14.40 ±0.12 


HD 91597 


Bllllne 




2000-02-04 


1.4 


2 


2030^= 














2000-02-06 


2.3 


2 


2100* 














2000-02-08 




4 


5629*= 









Note. — (a) At: time between successive observations. (6) Order of the polynomial used to fit the local 
continuum adjacent to the interstellar O VI feature, (c) Errors are la. For the equivalent widths, the first error 
is statistical, while the second one reflects the error in continuum placement. Wc did not allow for blending with 
the HD (6 - 0) i?(0) line at 1031.93 A. (d) from Prinja & Crowther (1998). (e) Only data from orbital night 
were included to reduce airglow contamination. (/) Voo from Howarth et al. (1997). 
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Rest, frame velocity (km s ') 

Fig. 1.— 

Top: FUSE spectra of Sk —69° 124 in the vicinity of O VI obtained at two times separated by 62 days. The rest wavelength and ai'c indicated 
for the O VI doublet by horizontal bars. The observed difi'erences between ~1027 and 1035 A are due to wind variability from both components of the 
doublet. The horizontal thick dotted line shows where the Galactic and LMC interstellar components would be expected. Bottom: Ratio of the fluxes 
plotted in the rest frame velocity of each component of the doublet. Note that while the wind variability does not affect the interstellar O VI A1038 
feature, the middle diagram shows clearly how the high velocity O VI A1038 wind variability can interact with the interstellar O VI A1032 feature. The 
dotted line indicates Vqo, while the extent of the variation is shown by the solid thick vertical lines. Error bars show the typical la error in the flux 
(top) or flux ratio (bottom). The positions of strong airglow lines due to H I and O I emission are indicated in both panels. Note that while the wind 
variability does not affect the interstellar O VI A1038 feature, the middle diagram shows clearly how the high velocity O VI A1038 wind variability can 
interact with the interstellar O VI A1032 feature. The dotted line indicates v^a, while the extent of the variation is shown by the solid thick vertical 
lines. Error bars show the typical Icr error in the flux (top) or flux ratio (bottom). The positions of strong airglow lines due to H I and o I emission 
are indicated in both panels. 
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Kcst frame velocity (km s"') Kcst frame velocity (km s"') liest frame velocity (km s"') 

Fig. 2. — Top and Middle: Continuum placement in the region of O VI A1032. Horizontal bars show the regions used to define the 
continuum. The dotted line corresponds to the continuum of the spectrum in the top panel, while the dashed line is the continuum of the 
spectrum in the middle panel. Bottom: Normalized fluxes resulting from the different continuum fits. The thicker line corresponds to the 
continuum from the middle panel. Dotted tick marks indicate the positions of H2 lines; solid tick marks show the positions of the O VI 
Galactic and LMC components (~ and 200 kms~^, respectively). The feature at ~ —135 kms~^ in the spectrum of HD 97913 is CI I 
A1031.3 
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Fig. 3.— Spectra of HD 91597 showing the propagation of a DAC in the O VI wind profile of A1032 (A < 1028 A) and A1038 (A < 1034 
A), corresponding to velocities of ~— 1170 kms~^ (day 2) and ~— 1370 kms~^ (day 3). The DAC accelerates slowly (~10~^ kms~^), which 
is typical of the motion of DACs observed in the winds of other stars (Prinja 1998). 



